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Abstract 

We show that the imphcit equation of a surface in 3-dimensional pro- 
jective space parametrized by bi-homogeneous polynomials of bi-degree 
{d,d), for a given integer d > 1, can be represented and computed from 
the linear syzygies of its parametrization if the base points are isolated 
and form locally a complete intersection. 

1 Introduction 

1 

Implicitization, i.e. finding the implicit equation of an algebraic curve or 
surface defined parametrically, is a classical problem and there are numerous 
approaches to its solution, most of them based either on resultants, Grobner 
bases, or syzygies. A good historical overview of methods based on resultants 
or Grobner bases can be found in [TSj and [Tj. 

Syzygy-based methods - also known as "moving surfaces" methods - were 
introduced in [13] and have been further developed in a number of publications. 
In the case of curves, these methods show that the linear syzygies of a given 
parametrization can be put together in a square matrix whose determinant 
is an implicit equation. Several attempts to generalize these results to the 
case of surfaces have been made. The construction of a square matrix whose 
determinant is the implicit equation requires the use of quadratic syzygies in 
addition to the linear syzygies. Moreover, these methods are only valid for 
birational parametrizations and need additional assumptions in the presence of 
base points. 
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Just to name some recent publications on syzygy-based methods, we cite [B], 
which treats the imphcitization of base-point-free homogeneous parametriza- 
tions and [1], which does the same for parametrizations with base points. In 
[T] a determinantal representation of the implicit equation of a bi-homogeneous 
parametrization is constructed with linear and quadratic relations, whereas |12) 
gives such a construction in the toric case. 

Recently it has been proved in [S] and [2| that surfaces parametrized by the 
projective plane can be represented and computed only using the linear syzygies 
of the parametrization, in the case where the base points are isolated and locally 
complete intersections. In some sense, this result is a natural generalization of 
the method of "moving curves" developed for planar curves in [13j : the only 
difference is that the matrix obtained in the case of surfaces is not square, but 
still represents the surface (see the end of Section [3] for a detailed explanation 
of this term). 

In this paper our main objective is to develop a similar imphcitization tech- 
nique for surfaces given by bi-homogeneous parametrizations, which are of in- 
terest for a number of applications in geometric modelling and computer-aided 
design. We will show that also in this case the surface can be represented by a 
non-square matrix constructed by only using linear syzygies and we will explain 
how to efficiently compute this matrix with standard computer algebra systems. 

More precisely, we focus on the following problem. Let K be any field (all 
the varieties we will consider hereafter are understood to be taken over K) . We 
suppose given a rational map 



where each polynomial /i, /2, /a, A is bi-homogeneous of bi-degree {d, d), d being 
a given positive integer, with respect to the homogeneous variables (s : u) and 
{t : v). We assume that 

• (j) parametrizes a surface (which is equivalent to require that is a 
generically finite map onto its image) which is hence irreducible 

• the greatest common divisor of /i, /2, /s, A is a non-zero constant which 
essentially requires the number of base points of (j) to be finite (possibly 
zero) . 

We aim to find a representation of in terms of linear syzygies of /i , /2, /s and 
fi similar to the known ones for plane curves and for space surfaces parametrized 
by the projective plane. 

The paper is organized as follows. In Section [5] we give an equivalent for- 
mulation of our problem which replaces the given N x N-graduation by a sin- 
gle N-graduation. In Section [3] we will introduce an associated approximation 
complex that will be used in Section 13.31 to prove our main result. Then an 
algorithmic version is detailed in Section 31 as well as an illustrative example. 



pi X pi 
(s : u) X {t : v) 
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2 The Segre embedding 



It is well-known that x P^ can be embedded in through the so-called Segre 
embedding 

pi X P^ A P3 

{s : u) X [t : v) ^ {st : sv : ut : uv). 

We denote by its image, which is an irreducible surface of degree 2 in P^, 
whose equation in the coordinates Xi,X2, X3, X4 of P^ is known to be X1X4 — 
X2X3. Our strategy to solve our implicitization problem is to reparametrize the 
surface by ^ C P^, that is to say to consider as the closed image of the 
map tp from ^ to ¥^ fitting in the commutative diagram 

pi X pi -^-^ p3 (1) 

P 

y 

In the rest of this paper we will use the map ip = (j) o p^^ to implicitize J^f, 
which has the advantage of replacing the N x N-graduation of P^ x pi by a single 
N-graduation. In order to justify this approach we need to describe explicitly 
the algebraic counterparts of the maps in the above diagram. 

We begin with the map 0. The polynomial ring ]K[s,zi] is canonically N- 
graded, 

IK[s, u] = ^ IK[s, u]„ = K[s, u]o © K[s, u\i © K[s, u]2 © • . • 

neN 

where IK[s,u]i denotes the degree i homogeneous component of IK[s,m], and its 
homogeneous spectrum is the projective line, i.e. Proj(K[s, u]) = Pjj. Of course, 
the same is true for the polynomial ring ]K[i,z;]. Now, consider the N-graded 
K-algebra 

5* := (K[s, u]n (g)K K[t, v]n) C K[s, u] (g,K K[t, v] 

neN 

which is finitely generated by Si as an iSg-algebra. Then P^ x pi is the homo- 
geneous spectrum Proj(S') of S. Introducing new indeterminates Ti, ^3, T4, 
the map is hence induced by the graded /c-algebra morphism 

K[Ti,T2,r3,T4] i S 

Ti fi{s,u,t,v) i = 1, . . . ,4. 

By Theorem 2.1], ker(p) C IK[Ti, Tj, T3, T4] is the defining ideal of the closed 
image of in P^ = Proj(K[Ti, . . . ,T4]); it is prime (since is a domain) and 
principal (since it is of codimension one by hypothesis and K[Ti, T2, T3, r4] is 
factorial), i.e. any generator of ker(p) gives an equation of . 
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We now turn to the Segre embedding p. As we did for the map (j) (note that 
the Segre embedding is itself a parametrization of a surface in projective space) 
the map p is induced by the graded fc-algebra morphism 

K[Xi,X2,X3,X4] ^ S 

Xi ^ st 
X2 ^ sv 
^ ut 
Xi ^ uv. 

However, in this case 9 is surjective and graded (it preserves the degree). More- 
over, it is easy to se£0 that its kernel is the principal ideal (X1X4 — X2X3) C 
'K[Xi,X2tX^,X4\. Therefore, 9 induces a graded isomorphism of N-graded K- 
algebras 

9:A:= K[Xi, X2, X^^, X^]/ {X^X^ ~ X^X^) ^ S 
which identifies x = Proj(5') with the Segre variety y = Proj(A) C P^ = 

Proj(K[Xi,X2,X3,X4]). 

We are now ready to describe ip. This map is of the form 

^ p3 ±^ p3 

(Xi : X2 : X3 : X4) ^ (51 : 52 : .93 : .94)(Xi, ^2, ^3, ^4) 

where 51,52,53,54 are homogeneous polynomials of the same degree in 
K[Xi,X2,X3,X4]. By the graded isomorphism 9, it follows that deg((/)) = 
deg(^) (we understand co-restriction to ^) and also that the 5i's must have 
degree d. To give an algorithmic construction we just have to determine the 
inverse map of 9. To do this, for all n G N define the integer /Cj-""* :— max(0, n — 
i — j) and consider the map 



s'u"-^i^i;"-J' ^ x^^'^''^^ X^^'"''-"' X^-"' 

(for all couples (i, j) £ {0, . . . , n}^). Then, we define the map 
:= u;„ : ^ -> K[Xi, X2, X3, X4] 



neN 



which induces the inverse of 9 by passing to A — IK[Xi, X2, ^3, -''^4]/(-'fi-'^4 ~ 
-'^2-'^3) (this is easy to check). Observe also that no monomial in the image of 



^ We clearly have {X1X4 — C ker(6). Now, if P G kcr(6') we deduce by a pseudo- 

euclidean division that there exists N N* such that 

Xf P = Q(Xi,. . .,X4){XiX4 - X2X3) + R{X2,X3,X4). 

But then R € ker(6) and it is obvious to check that we have K[X2, Xs, X4] n ker(6') = 0. 
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Lo is divisible by X1X4, so our representation of the inverse of B can be thought 
of as already reduced. Moreover, the coefficients of the /i's and the g^'s are in 
correspondence: only the monomials are changed by lo. 
Therefore, we proved 

Proposition 1 Defining for all i = 1, 2, 3, 4 the homogeneous polynomial 

g,{XuX2,X3,Xi) ■.^Lo{f,{s,u,t,v))eK[Xi,X2,X3,XiU, 

the map Q is a parametrization of the surface C with the property that 
deg(V') = deg((/>). 

Furthermore, we actually proved that our initial problem, namely the im- 
plicitization of <j) in terms of syzygies, is equivalent to the same problem with 
the parametrization -0 which is induced by the map 

K[Ti,r2,T3,r4] ^ A 

This can be summarized by the following commutative diagram, which is the 
algebraic translation of the diagram U) . 

e 

A 

This shows that the syzygies of the /^'s over S are in correspondence with the 
syzygies of the giS over A, in particular ker{h) = ker{p). Moreover, it also shows 
that the base points of the parametrization are in one-to-one correspondence 
with the base points of the parametrization -0 and that their local structure 
(complete intersection, multiplicity, etc.) is preserved by this correspondence. 

Another interesting remark is the following: By [Sj Theorem 2.5], we deduce 
that we have the equality 

deg(^)deg(jr) = deg{y)dP - ^ 

pev(ffi,...,gi)nj5^cp3 

where Cp denotes the algebraic multiplicity (in the sense of Hilbert-Samuel). 
Since it is immediate to check that deg(^) = 2 we recover the well-known 
formula of intersection theory (see Prop. 4.4] or [7l Appendix]): 

deg(0)deg(J^) = 2d2 _ ep. (3) 

pGy(/i,...,/4)cpixpi 

Therefore, in the rest of this paper we will focus on the implicitization of ip 
by means of linear syzygies, which is a completely equivalent problem to our 
initial one. 
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3 The approximation complex 



For simplicity, we will denote by Xi the classes of each variable in the quotient 
ringyl = K[2C\/ (X1X4—X2XS), where X stands for the sequence Xi, X2, X^, X4. 
Recall that A is canonically graded, each variable having weight 1. Let / = 
(ffii 52, 53, 94) C A be the ideal generated by the gi's. We give a brief definition 
of the approximation complex of cycles associated to the sequence gi, §2, 93, 94 
over A. This has been studied in depth in [IT], see also [14]. Under certain 
conditions this complex is a free resolution of the symmetric algebra Sym^(/), 
which is one of the main motivations for its study. Another essential feature of 
this complex is that - unlike the Koszul complex - its homology depends only 
on the ideal {gi, . . . ,174), not on the generators g^. Here is the construction: 

We consider the Koszul complex {K,{g, A),d,) associated to gi, . . . , 174 over 
A and denote Zi = ker((ii) and Bi = im(di+i). It is of the form 

A{-4d) ^ ^(-3d)4 ^ Ai-2df ^ Ai-d)^ ^ A 

where the differentials are matrices whose non-zero entries are ±171, . . . , ±174. We 
introduce new variables Ti, . . . , T4 and set Zi = Zi{i ■ d) ^[2^], which we will 
consider as bi-graded ^[T]-modules (one grading is induced by the grading of A, 
the other one comes from setting deg(Ti) = 1 for all i). Now the approximation 
complex of cycles (Z, (g, A), e,), or simply Z,, is the complex 

- Z3(-3) ^ Z2{~2) X Zi(-l) ^ Zo 

where the differentials e, are obtained by replacing gi by Ti for all i in the 
matrices of d, (note that Z4 — 0, since d4 is injective). It is an important 
remark that 

im(ei) = Iy^P^T^^ I e A[T],Y,P^9^ = o| (4) 
1 1=1 1=1 J 

= \ T.P^^- I ^ A,Y,p,g, = 0) C A[T] 
\i=i i=i / 

and therefore Ho{Z,) = A[T]/im(ei) ~ Sym^(J). Note that the degree shifts 
indicated in the complex above are with respect to the grading given by the T^'s, 
while the degree shifts with respect to the grading of A are already contained 
in our definition of the Zj's. From now on, when we take the degree v part of 
the approximation complex, denoted {Z,)^, it should always be understood to 
be taken with respect to the grading induced by A. Hereafter we denote by m 
the ideal (Xi, X2, X3, X4) C A. 

3.1 Acyclicity criterion 

Our first concern is to show that Z, (171, . . . , 34; A), the approximation complex 
of cycles is acyclic under suitable assumptions. We have, similarly to |3, Lemma 
2], the following 
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Lemma 1 Suppose that I = [gi, 92, 93, §4) C A is of codimension at least 2, 
and let ^ Proj(A//) C -y. Then the following are equivalent: 

(i) Z, is acyclic, 

(a) is acyclic outside V^(tn), 

(Hi) ^ is locally defined by 3 equations (i.e. locally an almost complete inter- 
section). 

Proof. The proof is very similar to [31 Lemma 2]; the only difference is that A 
is not a polynomial ring here, but it is still a Gorenstein ring which is the main 
required property for A. Observe that the lemma is unaffected by an extension 
of the base field, so one may assume that K is infinite. 

By [m Theorem 12.9], we know that Z, is acyclic (resp. acyclic outside 
y(m)) if and only if / is generated by a proper sequence (resp. ^ is locally 
defined by a proper sequence). Recall that a sequence ai, . . . , a„ of elements in 
a commutative ring B is a proper sequence if 

ai+iHj{ai, . . . , tti; B) — for i = 0, . . . , n — 1 and j > 0, 

where the Hj 's denote the homology groups of the corresponding Koszul com- 
plex. 

It is clear that (i) implies (ii). Assuming (ii), we will now deduce that ^ 
is locally defined by a proper sequence. As explained in 3", Lemma 2], one can 
choose /ii, ft,2, /13, /i4 to be sufficiently generic linear combinations of the gt's such 
that 

• {hi, . . . ,/i4) = (gi, ... ,34) C A, 

• /ii,/i2 is an A-regular sequence, which implies that hi, ft,2, is a proper 
sequence in A, 

• hi, . . . ,hi form a proper sequence outside V{m). 
By [21 Theorem 1.6.16], we have 

Hi{hi,h2,h3;A) ~ Ext\ {A/{hi,h2,hi), A) 

and since A is Gorenstein (for it is a complete intersection), i.e. isomorphic to 
its canonical module [2 Theorem 3.3.7], then 

i/i(/ii,/i2,/i3;^) -Ext^(A/J,^) -w^/j (5) 

outside V{m), where cij_ stands for the canonical module and J :— {hi, /i2, /is) C 
A. Since the annihilator of oja/j over A is {J : m°°) C A (observe that A/ J 
defines isolated points and use for instance [SI Corollary 21.3]), we deduce that 
ft.4 G ( J : tn°°), that is to say that 3^ is locally defined by 3 equations. 

Now, assume (iii). Similarly to what we did above, one can find hi, . . . ,hi 
sufficiently generic linear combinations of the giS so that hi, /i2 is an ^-regular 
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sequence and /ii, ft.2, ^3 define It follows that /14 G ( J : m°°) C A, where J := 
(/ii, /i2, C A, and hence (O implies that ft.4 annihilates Hi{hi, h2, h^] A); it 
follows that /ii, . . . , /i4 form a proper sequence in A, so 2^, is acyclic. □ 

As soon as the base points (if there are any) of the parametrization ip 
(or equivalently (j)) are isolated and locally defined by 3 equations, then its 
associated approximation complex of cycles is acyclic. Therefore, it can be 
used to compute and represent the codimension one part of the annihilator of 
the A[Ti,. . . ,T4] -module H^{Z,) which is nothing but the symmetric algebra 
Sym^(/). Our interest in this module is motivated by the following 

Lemma 2 Suppose that £^ :— Proj(yl//) has dimension < and is locally 
defined by 3 equations. If rj is an integer such that 

i/^ (Sym^ (/)), = foralli^>rj, 

then, for all v >rj we have 

annK[T](Sym^ (/),,) = anuKf^] (Sym^(/)^) C ker{h). 

Moreover, the above inclusion is an equality if ^ is locally defined by 2 equa- 
tions. 

Proof. For all > ry, the equality 

annK[7](Sym^(/)i.) = annK[T](Sym^ (/),,) 

is proven in [H Proposition 5.1] for A — K[X]. However, the same proof can 
be applied without modifications to our setting: The key property used in the 
proof is the fact that the canonical map Ai ® An — > An+i is surjective and this 
is also vahd for A = K[X]/{XiXi - X2X^). Moreover, by (g]) we have that 
annK[T](Sym^(/)^) / for ^ if and only if ^ is locally generated by at 
most 3 equations, and in this case it is clear that it is contained in ker(/i). Fi- 
nally, if 3^ is locally defined by at most 2 equations, meaning that 3^ is locally 
a complete intersection, then / is of linear type outside V{m) (use for instance 
[51 Propositions 4.1 and 4.5]) which shows the last claimed equality as proven 
in O Proposition 5.1]. □ 

In other words, if the base points of the parametrization are isolated and 
locally complete intersections then certain graded parts of the approximation 
complex Zt yield a way to compute an implicit equation of J^. Our next task 
is to explicitly describe the saturation index of the symmetric algebra, i.e. the 
integer rj appearing in Lemma [2] This will provide us with the key tool for 
developing the algorithm presented in Section HI 

3.2 The saturation index 

For any ideal J of yl we denote by J^^' the saturation of J with respect to the 
ideal m, i.e. J^^* := (J :a m°°) C A. Also, we recall that if M is a N-graded 
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S-module, where _B is a N-graded ring, its initial degree is defined as 

indeg(M) := min{i^ £ N : il/^ 7^ 0} > 0. 
With these notations, we have 

Theorem 1 Ij :~ Proj(A//) is a zero- dimensional scheme (i.e. supported 
on a finite number of points, possibly zero) then 

i/,'i(Sym^(/)), = Vi. > 2d - 1 - indeg(r'^*). 

The proof of this theorem is actuaUy similar to the proof of [3l Theorem 
4]. The difference is that in our case the ring A is not a polynomial ring but 
a quotient ring. So to validate the proof of f3, Theorem 4] we have to make 
explicit the local cohomology and the dualizing module of A which is, as a 
complete intersection, a Gorenstein ring (the key property for what follows). 
We state these results in a little more general case for the sake of clarity. 

Proposition 2 Let k be a commutative Noetherian ring and C := k[Xi, . . . , Xn], 
with n > 1, which is canonically graded by deg(Xi) — 1 for all i = 1, . . . , n. Sup- 
pose given a homogeneous polynomial f of degree r > 1 and consider the graded 
quotient ring B := C/{f). The following properties hold: 

• ujb — B(—n-\-r), a graded isomorphism where ujb stands for the canonical 
module of B, 

• (B) — if i ^ n — \ and for all E Z 

Hl~\B),^B{-n + r).,, 

• if denotes the Koszul complex associated to a given sequence (ai, . . . , as) 
of homogeneous elements in B of degree di,...,ds respectively, then we 
have the isomorphisms 

s 

Hl~\K,), ~ HomB/„(/C-.(^d» -n + r)_„B/m). 

i=l 

Proof. To prove the first claim, we first recall that we have loc — C{~n). Then, 
[H Corollary 3.6.14] shows that 

UJB ^ {loc / f ■^c){r) - B{-n + r). 

For the second claim, we recall that the local cohomology of C is well-known: 
F4(C) = for allien and 

Hl{C),^C-n-. (6) 
for all V eZi. Now, the exact sequence 
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whose long exact cohomology sequence contains the segments 
Hi{C) Hl{B) Hi+\C{~r)) 



imphes that (B) = for aU j < n — 1 as for j + 1 < n both the left and the 
right hand side vanish. Furthermore, the segment 

taken in degree f shows 

i/r' (B)^ = ker[Hl{C{~r)), Hl{C),) . 

By the self-duality of the Koszul complex and ([6]) this later equals exactly 
B-u-n+r- Finally, since dim(C) = n we have dim(i?) = n — 1 which implies 
that Hi{B) = for j > n - 1 by [2, Theorem 3.5.7]. 

The third claim is a direct generalization of the classical property 

s 

- n)_^, C/m). 

The only thing which changes is the shift by r in the canonical module of B and 
the dimension of B which is n — 1 whereas dim(C) = n. □ 

Proof of theorem [TJ We consider the two spectral sequences associated 
to the double complex H^{Z,): 



CO (Z3) ^ Cl{Z^) C-o (Zi) ^ Cl{Z,) 

^ Cl{Z^) Cl{Z2) Cl{Z,) Ci(Zo) ^ 
\ \ \ \ 

^ C4 (Z3) CiiZ^) (Zi) ^ Ci{Z,) 

They both converge to the hypercohomology of Z,. One of them stabilizes at 
level two with: 

r Hl{Hq{Z,)) forp = 0,l andg>0 
2'^^ = J El - <^ i/S(Sym^(/)) for q = 
[ else. 

and the other one gives at level one: 

,"El = Hl{Z,)[qd]®AA[T]{~q). 
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As explained in [31 Theorem 4] , the comparison of these two spectral sequences 
and [21 Lemma 1] sho'wH that the module i/^(Sym^(/))i, vanishes as soon as 
{i"EP)i, vanishes for p — 2,3. Moreover, setting — * :— Homgr^( — ,A/m), we 
have the graded isomorphisms 

i"El^{A/iy[2~d\^AA[T]{-3) 

and 

^ (r^V/)*[2 - 2d] d) A[T]{-2). 

It follows that [i" £2)1, and {i" E^)^ vanish simultaneously if 

V > min(d - 2, 2d - 2 - indeg(/"''V/)). 

This is true whenever v >vq :— 2ci— 1— indeg(/''^*), since min((i, indeg(/''^*//)) = 
indeg(/''^*). □ 

Remark 1 Since I is generated in degree d and I C P'^^ we have the inequality 
< indeg(r'»^*) < d. It follows that 

d - 1 < 2rf - 1 - indeg(r^*) < 2d - 1. 

The lower bound is reached whenever the ideal I is saturated (meaning I = 
P^^J and the higher bound corresponds to the abscence of base points of the 
parametrization. 

3.3 The main result 

We now have all the tools necessary at our disposal and can proceed to the 
main result of this paper. But before, recall that there are two distinct notions 
of multiplicity for a base point p G V{I) C\ .y Cf'^: the algebraic multiplicity 
denoted Cp and the geometric multiplicity denoted dp (see for instance §2.2] 
for more details). 

Theorem 2 Assume that dim ^ :— Proj(A//) < and that 3^ is locally an 
almost complete intersection (i.e. locally defined by 3 equations). Then, for 
every integer 

v>uo:=2d-l- indeg(/"^') 

the determinant D of the complex {Z,),j of K.\T\-modules (which is unique up to 
multiplication by a non-zero constant in M^) is a non-zero homogeneous element 
in K[T], independent of v > vq and of degree 

2d^ - J2 '^P 

p£V{i)n,y(iP'-> 

such that D = F'^'^'^^^'^G where F is the implicit equation of M' , G is coprime 
with F and dcg(G) = Y.peV{i)nA^P - dp). 

Moreover, G € K \ {0} if and only if is locally a complete intersection 
(i.e. locally defined by 2 equations). 

^Note that [3] Lemma 1] can be applied verbatim in our case (modulo some little change 
on the degree shifts that we will describe below) because of Proposition [2] 
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Proof. First of all, observe that D is independent of v by theorem [TJ It is 
an homogeneous element of K[T] because (2.)^ is a graded complex of K[T]- 
modules and it is non-zero because £P is locally an almost complete intersection, 
a fact we already used in Lemma [5] 

The computation of deg{D) can be done as in pT, Theorem 4]: For ^ we 
have 

deg(L») = dim(Zi)^+d - 2 dim(Z2),y+2d + 3 dini(Z3),,+3d. 

In the case where all the i?i's, with i > 0, vanish then deg{D) — 2<P . If Hi and 
H2 are non-zero, then they contribute to the above quantity for 

dim(i/i),y+d - dim(iJ2)i/+d - 2dmi{H2)u+2d 

= dim(i7o)i.+d - 2 dim(iJ2),.+2d = -deg^ (7) 

where we assume that v ^ since H2 — i-OAji- Therefore, we deduce that 

degCL*) = 2d^ - deg^ = 2d^ ~ ^ dp. (8) 

p(EV(/)nj^cP3 

Now, setting q :— ker(/i) and using standard properties of determinants of 
complexes we compute 

[dct((Z.).)] = div(iJo(2.)) 

= div(Sym^(/),,) 

^ length((Sym^(/).)p) • [p] 

p prime, codim(p) — 1 

= length((Sym^(/),),)-[q] + -- - . 

Since length((Sym^(/)^)q) = deg(V') as proved in [5] Theorem 5.2], we deduce 
that D = F'i<=s(''')G where G does not divide F. 

Finally, using equations ([3]) and ([8]) we deduce that 

deg(G)= (ep-dp), 
pev(/)nj5' 

and it is well-known that Cp > dp with equality if and only if the point p is 
locally a complete intersection. □ 

Recall that the determinant of the complex (2,)^ can either be obtained as 
an alternating product over some sub-determinants of the matrices appearing in 
the complex or as a gcd of maximal minors of the first map in the (2,)y-complex 
(we will explicitly construct this matrix M in the next section). One can either 
compute this gcd directly or factorize one of the maximal minors, however, both 
methods are computationally expensive (as all existing implicitization methods). 
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From a practical point of view, it might be interesting to avoid the actual 
computation of an implicit equation and use instead the matrix M as an im- 
plicit representation of the surface, since it is more compact and much easier to 
compute. To give an example, let us suppose that we are in the case of locally 
complete intersection base points. Then if we want to decide if a given point P 
lies on the surface there is no need to compute the implicit equation: It suffices 
to evaluate M in this point, as the rank of M drops if and only if P belongs to 
the surface. 

This is due to the fact that for a commutative ring R and a morphism 
a : i?" with m > rt we always have 

ann7j(coker(Q;))" C /„(q;) C ann^f (coker(Q!)) 

where denotes the ideal generated by the maximal minors of the matrix 

of a, i.e. the principal Fitting ideal of a (see for instance [U Proposition 20.7]). 
Ours is the special case R — K[T] and a is the first map in i.e. the one 

induced by ei, and hence coker{a) — Sym^(/),y. Geometrically, this means that 
the maximal minors of M define the hypersurface Jff by Lemma [21 and conse- 
quently, the points for which the rank of M drops are exactly those belonging 
to Jf. 

Similarly, other problems arising from applications might be solved by direct 
computations using the matrix representation without the (expensive) transition 
to the implicit equation. 

4 Algorithm 

In order to show explicitly how the theoretical results from the previous sections 
are used in practice, we formulate an algorithm for the actual computation of the 
matrix representing the implicit equation. It is efficient and easy to implement, 
as it consists basically of the resolution of a linear system. We give only the 
essential steps, see (Sj Section 3] for a more detailed description of a very similar 
algorithm. 

• Given four bi- homogeneous polynomials fi, f2, fs, fi of degree d, define 
the homogeneous polynomials 51,32, 53,54 E A = K[2C\/{XiX4 ~ X2X3) 
of the same degree by setting 5^ = oj{fi), where lo is the isomorphism 
defined in Section [21 

• Find the solution space W of the linear system (over K) defined by 

-*e{i,...,4} 

where (01,02,03,04) G {A^^^)^ and = 2d — 1 — indeg(P''*), i.e. one 
writes the equation with respect to a basis of A^^^d and compares the 
coefficients. W is represented by a diia^X^va+d) x 4dimK(Aj,Q)-matrix TV, 
where the first k :— Ami-^{Ay^^) columns represent the coefficients of oi, 
the next k coefficients 02, etc. 



13 



• For i e {1, . . . , 4}, let Afj be the k x /c-matrix Ti ■ Idfe. Then 



/Ml 



M :=N- 



is a matrix of the first map of the graded part (Zt)^^,of the approximation 
complex. 

As we proved, in the case where the base points of the parametrization (f) are 
isolated and locally complete intersections, M represents the surface J^. Also, 
the gcd of the maximal minors (of size k) of M equals its implicit equation. 



An illustrative example 

We now present an example to illustrate our method, which provides a matrix- 
based representation of the implicit equation of by means of the linear syzy- 
gies of its parametrization (j) (or more precisely, of ip)- It should be emphasized 
that all the following computations are presented in order to explore in detail 
our approach and are not all required to get the expected matrix-based repre- 
sentation. Our code is written for Macaulay2 (see [10]), in which one can easily 
compute all the terms and maps of the approximation complex. 
Consider the following example taken from [1] Example 4.16]: 

S=QQ[s,u,t,v] ; 
d=2; 

f l=u~2*t*v+s~2*t*v 
f 2=u~2*t~2+s*u*v"2 
f3=s~2*v~2+s~2*t"2 
f4=s~2*t*v 

F=matrix{{f 1 , f 2 , f 3 , f 4» 

Note that the interested reader can experiment with his own example just by 
changing the above definitions of the polynomials /i, /2, /s, /4 giving the para- 
metrization. 

The first thing to do is to use the isomorphism 6 to switch from S (note 
that the ring S defined in the above command is not exactly the ring S we have 
introduced in Section [2]) to the ring A: 

SX=S [xl,x2,x3,x4] 
F=sub(F,SX) 
ST={}; X={}; 
for i from to d do ( 
for j from to d do ( 
k=max (0 , d-i- j ) ; 

ST=append(ST,s~i*u" (d-i)*t~j*v~(d-j)) ; 
X=append(X, 
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xl- (i+j -d+k) *x2- (d-j -k) *x3- (d-i-k) *x4- (k) ) ; 

) 

) 

ST=matrix {ST}; ST=sub(ST,SX) ; X=matrix {X}; 
(M , C) =coef f icients (F , Var iables=> 

{s_SX , u_SX , t_SX , v_SX} , Monomials=>ST) 
G=X*C — this is the parcunetrization, but in SX 
A=QQ [xl , x2 , x3 , x4] / (xl*x4-x2*x3) 
r=map (A , SX , {xl , x2 , x3 , x4 , , , , 0}) 
G=r(G); 

G=matrix{{G_(0,0) ,G_(0,1) ,G_(0,2) ,G_(0,3)}} 

The matrix G is the matrix (with entries in A) of the parametrization ^ from the 
Segre variety S^tof^. One should note that the quotient ring yl is a very simple 
quotient ring: essentially, computations in A can be done in K[Xi,...,X4] 
modulo the substitution of by Moreover, bases for A in any 

given degrees can easily be pre-computed since they do not depend on the given 
parametrization (p. 

We can now define the terms of the approximation complex of cycles : 

Z0=A-1; 

Zl=kernel koszul(l,G); 
Z2=keriiel koszul(2,G); 
Z3=kernel koszul(3,G); 

As we already remarked, Z4 = 0. Define the integer 
nu=2*d-l 

We can compute the Euler characteristic of (Z,)^ and check that it is zero with 
the command 

hilbertFunct ion (nu , ZO) -hilbertPunction (nu+d , Zl) + 
hilbertFunction(nu+2*d, Z2)-hilbertFunction(nu+3*d, Z3) 

and also compute the degree of D, the determinant of the complex {Z,)iy, with 
the command 

hilbertFunction(nu+d, Zl) -2*hilbertFunction(nu+2*d, Z2) 
+3*hilbertFunction(nu+3*d, Z3) 

This number equals the degree of if all the base points, if any, form locally 
a complete intersection. In this example, we find degree 7. 

At this step, one can try to lower the integer i/ according to Theorem [2l to 
this end we compute the degrees of the generators of the saturation of the ideal 
{91, ■■■,94)- 

degrees gens saturate (ideal G, ideal (xl,x2,x3,x4)) 
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Since, in this example, the smallest degree is 1 we can redefine 
nu=2*d-2 

and we can re-check the Euler characteristic and the degree of the determinant 
of (Z.),. 

We can now compute the matrix of the first map of {2,)i,, that is to say the 
matrix of linear syzygies of ffi, • • • , 54 which represents Jf : 

R=A [T1,T2,T3,T4] 
G=sub(G,R) ; 

Zlnu=super basis(nu+d,Zl) ; 

Tnu=inatrix{{Tl , T2 , T3 , T4}}*subst itute (Zlnu , R) ; 
(in,M) = 

coef f icients (Tnu , Var iables=>{xl_R , x2_R , x3_R , x4_R} , 
Monoinials=>substitute(basis(nu,A) ,R)) ; 

The matrix M is the desired matrix, and it is of size 9 x 12. 

5 Comments and conclusion 

We have presented a new approach to compute an implicit representation in 
terms of linear syzygies for a surface in P'^ parametrized by bi-homogeneous 
polynomials of bi-degree {d,d), d > 1, under the assumption that the base 
points are isolated and locally complete intersections. This result, along with 
the similar ones for parametrizations over the projective plane, shows that in 
many cases it is not necessary to use quadratic syzygies in order to represent 
the implicit equation of a surface. 

We should point out that this method has the advantages of being valid in 
a very general setting (we have neither assumed birationality nor made other 
additional assumptions on the paramctrization) and of working well in the pres- 
ence of base points. Furthermore, the matrix representing the surface can be 
computed in a very efficient way. 

It would be nice if we could use the same method for mixed degrees as 
well, i.e. consider parametrizations by bi-homogeneous polynomials of bi-degree 
(di, ^2) with di,d2 > 1. Let us discuss some ideas on how to generalize to the 

mixed case: 

• Putting weights on the variables in S will not give us good properties for 
S, for instance S will not be generated by Si as an So algebra in general. 

• Considering the bi-degree (max((ii, ^2), max(di, c?2)) is not possible be- 
cause it introduces a base point locus of positive dimension and we will 
lose the acyclicity of the approximation complex. 

• One way to come back to unmixed bi-degree is to make the substitutions 

g ^ glcni{di,d2)/di g^j^^ ^ ^ ^Icm(di,d2)/<i2 
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Everything works fine in this case, but we are not representing F'^'^^^'^\ 
but Fdog(v)icm(di,d2)/gcd(di,d2) ^^[(.-[^ is is not Optimal, as it increases the 
size of the matrices involved. For instance, we could treat Example 10 
from [12] in this way. It is a surface of bi-degree (2,3) defined by 

/i = {t + t^){.s-lf + {l + .st-sh){t-lf 

f2 = i-t - t^){s - If + {-I + St + Sh){t - 1)^ 

h = {t-t^){s-lf + [-l-st + sH){t-lf 

/4 = {t + e){s-if + St- sh){t-if 

By replacing s by s^ and i by t^, we obtain a parametrization of bi-degree 
(6,6) and can indeed be computed in degree i/>2-6 — 1 — 6 = 5of the 
approximation complex as the gcd of the maximal minors of a 42 x 36- 
matrix, whereas in the original paper it was computed as the determinant 
of a 5 X 5-matrix. 

Therefore, it seems that the tools we used above (and which work well for 
unmixed bi-degree) are not well-suited for this more general case and that it 
might be necessary to take into account the bi-graded structure of S in order to 
devise a method that is adapted to mixed bi-degrees. We hope to develop this 
in the near future. 
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